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Abstract 

We present a microscopic calculation of transport in strongly doped super- 
lattices where domain formation is likely to occur. Our theoretical method 
is based on a current formula involving the spectral functions of the system, 
and thus allows, in principle, a systematic investigation of various interaction 
mechanisms. Taking into account impurity scattering and optical phonons we 
obtain a good quantitative agreement with existing experimental data from 
Helgesen and Finstad (J. Appl. Phys. 69, 2689, (1991)). Furthermore the 
calculated spectral functions indicate a significant increase of the average in- 
tersubband spacing compared to the bare level differences which might explain 
the experimental trend. 

I. INTRODUCTION 

In semiconductor superlattices the electric transport is dominated by resonances between 
the localized energy levels inside the wells. A resonance structure appears in the current- 
voltage characteristics if the energy levels of different wells align Between these res- 
onances regions with negative differential conductivity (NDC) are likely to appear. This 
yields interesting current- voltage characteristics as already shown by Esaki and Chang ||. 
If the voltage applied along the growth direction of the superlattice yields an average field 
which is in the NDC region, the homogeneous field distribution breaks up and electric field 
domains form causing many current branches in the current-voltage characteristic (almost 
equal to the number of quantum wells). This has been extensively studied experimentally 
during the last decade 

Theoretically the complicated measured I(U) characteristics could be reproduced using 
models based on the combination of rate equations between the wells and Poisson's equation 
TQ-Fj . These models used simplified expressions for the tunnelling current between the 



wells |nj, or a fitted local current-field characteristics [11|. A theoretical approach for 
the transport in superlattices with strong coupling between the wells has been reported 



in Ref. |L3| with the restriction to one miniband. This restriction is not appropriate for 
situations where high field domains occur, as they are close to the resonance between the first 
and the second subband ||. Nevertheless, no quantitative calculations have been performed 
so far. 



The aforementioned theories for domain formation showed that the electric field as well as 
the carrier density are almost constant within the field domains. Thus the current density is 
determined by intrinsic features of the superlattice. This is a much simpler situation than in 
many other semiconductor structures, such as the double-barrier resonant-tunnelling diode, 
where the contacts strongly influence the field profile and the current densities by carrier 
injection. Furthermore the single branches extend up to a current density, where the low field 
domain reaches the maximum of the current-field relation. Therefore the current maxima of 
the branches are a good estimate for the maximum in the current-field relation, which does 
not depend strongly on the contacts. The situation becomes slightly more complicated, if 
the superlattice exhibits fluctuation in its parameters, which have a strong impact impact 
on the individual branches fl4] . 

In this paper we develop a unified microscopic theory, without adjustable parameters, 
to describe the underlying current-field relation. In contrast to the situation considered 
in Ref. flSfl , we consider the situation of weakly coupled quantum wells (Hi <C T) where 
Hi is the coupling between the wells and T is the broadening of the levels within the single 
wells. This means that the intrawell scattering rate is larger than the tunnelling rate. This is 
relevant to most experiments including that described in Ref. with which we will compare 
our results. 



II. THE MODEL 

We consider weakly coupled semiconductor quantum wells of period d. Then the electrons 
are localized in the wells and a reasonable basis set of wave functions is given by a product of 
Wannier functions ^> v (z — nd) localized in well n, and plane waves e 1 --. Here the z direction 
is defined to be the growth direction and k, r are vectors within the (x, y) plane, v denotes 
the subband within the well. 

Restricting ourselves to the lowest two minibands (denoted by a and b) and coupling 
between neighbouring wells we consider the following hamiltonian (F is the electric field, 
and e < is the charge of the electron) : 



H = £ [E»(k)al (k)a n (k) + E b n (k)bi(k)b n (k)] (1) 

n,k 

Hi = J2 [Ti a ai +1 (k)a n (k) + T b b+ + i(k)b n (k) 

n,k 

-eFRfai +1 (k)b n (k) - eFR ba bl + i(k)a n (k) + h.c] (2) 
H 2 = Y, [-eF(R a b ai(k)b n (k) + R b %(k)a n (k))} (3) 

n,k 

with E^(k) = E u +h 2 k 2 / (2m w )—eFn (m w is the effective mass in the well), and the couplings 
B£> = $ dzW' (z-hd)zW(z). 4|Tf| and A\T\\ are the miniband widths of subband a and 6, 
respectively. The term H 2 can be incorporated into the one electron states by diagonalizing 
ifo + H 2 [0. This leads to renormalized coefficients in Hq and Hi but does not change the 
structure of the problem for a homogeneous electric field. In a similar way higher subbands 
and continuum states will change the coefficients, so we prefer to omit H 2 in the following. 
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The total current from subband v in well n to subband \i in well n + 1 can be described 
by the following expression which is derived in Sect. 9.3 of Ref. [15|: 



v-> v 2 f°° d£ , 



,, h (n+l)k,nk'\ J_ 00 2llh 

■A% +1 (k, £ + fi n - fi n+ i) [n F {£) - n F (£ + fi n - Hn+i)} . 

Here \x n is the electro-chemical potential in well n and n F {£) = (1 + e 136 )" 1 is the Fermi 
function. The energy £ is measured with respect to \i n . Note that for equal densities in 
both quantum wells we have \x n — fi n +i — eFd. 

An{k,£) denotes the spectral function for the state k of the subband v in well number 
n. It is calculated in equilibrium neglecting the coupling to the other wells and is related to 
the retarded self-energy S^ ret (fc, £) via 

— 2ImS vret 

A n (t, £) - — - — - e ^ - ReS , rct) 2 - ( ImS ,rct)2 ( 5 ) 

If no perturbation is present the state k has a fixed energy, E^(k), and the spectral function 
becomes a ^-function A v n {k, £) = 2n5(£ + /x n — E^(k)). If scattering is present the states k 
are no longer eigenstates of the full hamiltonian and the spectral function becomes smeared 
out, which is often modeled by a Lorentzian using a constant S^ ret . 

While the full derivation is slightly tedious [|HJ the formula (f|) can be motivated 
quite easily: In the long-time limit energy has to be conserved during transitions caused 
by the time-independent interwell couplings H^ +l , knk ,. Therefore we have to consider 
tunnelling processes for a certain energy £ and integrate over £ afterwards. The factor 
[np{£) — np{£ + \i n — /i n +i)] takes into account the thermal occupation at the given energy 
in both wells. The free particle state &/ has a weight A n (ti , £) / (2it) in well n. Its transition 
probability to the state k in well n + 1 is given by 2-K\H l ^' +l ^ k nk ,\ 2 / h. The final state has a 
weight An + i(k, £ + \i n — yU n+1 )/(27r) at the given energy. Obviously one has to sum over all 
free particle states k, fc'. Finally, the factor 2 is due to the spin degeneracy. 

We calculate the self-energy £ aret (7c, £) for the lower subband for impurity scattering 

against the ionized donors within the self-consistent single-site-approximation (SSA) (which 

contains all noncrossing diagrams as shown in the inset of Fig. |l]). This approximation 

was also used in Ref. ]16|] to calculate the spectral functions in a single quantum well. 

The matrix element Vh(p) of the interaction is calculated from the given Wannier functions 

both for impurities in the same well (h = 0) and for remote impurities located in the h- 

th well counted from the electron, which become important for pdh < 1. We consider 

screening by the free 2D electron gas located both in the same well and in different wells 

within the T = Random Phase Approximation. The screened matrix element for impurity 

scattering is denoted by V£ c (p). Within the Born approximation the total scattering rate 

— 2 
with momentum transfer p is proportional to J2h ^/TG°) • m Fig. we have plotted this 

quant i t y as well as Vq c (p) as a function of p. For comparison we have also shown the 
result for the case when screening and scattering is restricted within the same well; this is 



larger than Vq c (p) because of the weaker screening. From Fig. |T] we can see that even for 
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p = the main contribution to the total scattering J2h 

2 



V£ c (p) comes from the scattering 



within the well 



v sc (p) 



Thus we can restrict ourselves to this contribution in our case. 



The self-consistent Born- Approximation (BA) (consisting of the first diagram shown in 
the inset of Fig. [[]) breaks down for matrix elements ircV ~ 1, where c = m w /(2nh 2 ) is the 
density of states per spin of the free 2D electron gas. Fig. [I] tells us that we are in the range 
where the higher order diagrams become important which invokes the need for the SSA. 

For the calculation of the self-energy E 6ret (fc, S) of the upper subband we additionally 
include the emission of optical phonons for scattering between the upper to the lower subband 
describing the intersubband relaxation. This leads to an additional imaginary part — Tp h in 
the self-energy E bret (fc, £). 



III. RESULTS OF THE CALCULATION 



Here we apply our approach to the experimental situation of Refs. |6|,|T7||. The 



Alo.3Gao.7As superlattice used there had a nominal barrier width b = 12 nm and a well 
width w = 8 nm. The middle 7 nm of the wells were n-doped with a doping density 
N D /A = 8.75 ■ 10 15 /m 2 . We calculate the coefficients in Eqs. (0,||,|D within the Kronig- 
Penney model assuming a parabolic dispersion with the conduction band offset AE C = 0.24 
eV, and the effectives masses m w = 0.067m e and my, = 0.0919m e for well and barrier, respec- 
tively |]I8" |. The resulting parameters are given in Table |I|. Just as in the more sophisticated 



calculation used in Ref. @JTj the calculated subband spacing E b — E a = 117.9 meV is smaller 
than the experimental value 124 meV determined by intersubband absorption. 

We approximated the doping profile by two 5-doping layers located at a distance of 
1.75 nm from the middle of the well. The calculated spectral functions for subband a are 
given in Fig. |2] for both the BA and the SSA. For Ej. = 0.02 eV both spectral functions 
A a (Ek, E s ) (where E s = £ + // — E a is measured with respect to E a ) are quite similar and 
have approximately a Lorentzian shape with a half width of T a « 0.01 eV. In contrast 
to this for Ek = the spectral function exhibit significant differences for the BA and the 
SSA. Within both approximations they exhibit a sharp onset at E s ~ —0.02 eV and are 
significantly different from Lorentzians. 

The broad width of the spectral functions has obviously an impact for the optical ab- 
sorption a(uj) between the subbands which is proportional to 

/oo 
dSn F (S)A a (k, E)A\k, £-hw). (6) 

Our calculation reveals a maximum at hu = 123.5 meV which is in excellent agreement with 
the experimental value. Nevertheless the full width at half maximum (35 meV) is larger 
than the experimental value 19 meV. (Within the BA the maximum is at hu = 123 meV 
and the width is 30 meV.) This shows that the measured interssubband spacing may deviate 
significantly form the bare energy levels for strongly doped samples. 

For the given density of carriers provided by the doping we have calculated the chemical 
potential for the actual density of states determined from the spectral functions. 

Finally we compute the current from Eq. (|4j). We use zero temperature and the Fermi 
functions n F (£) become step functions 0(— £). Then only the lowest subband is occupied 
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in equilibrium and the total current is given by the sum I = J a ~* a -f I a ^ b . Using the 
nominal sample parameters we obtain current maxima for the current-field characteristics 
which are almost an order of magnitude lower than the currents measured in the respective 
regions. This indicates that the coupling between the wells should be stronger. Using the 
barrier width b = 10.5 nm we find good agreement with the experimental data. In their 
own theoretical analysis of the data Helgesen, Finstad, and Johannessen [|17| use the width 
b = 10.8 nm, which indicates that they consider this value to be within the experimental 
uncertainty. The result of the calculation for b = 10.5 is given in Fig. [| The I(eFd) relation 
exhibits a first maximum at eF^jd = 13 meV with a maximum current of Im = 0.592 mA 
and a minimum at eF m d = 69 meV with a current I m = 0.086 mA. At eFd = 123 meV 
there is a second maximum due to the resonance between the first and the second subband 
with a maximum current / = 29.4 mA. The presence of a range with negative differential 
conductivity for 13 meV< eFd <69 meV causes an instability leading to the formation of 
field domains which is discussed in detail in Ref. [I9~| for an arbitrary v (F) relation. The 



domain branches in the characteristic should exhibit a maximum current Im if the high field 
domain is located at the anode and a lower current if the high field domain is located at the 
cathode, which concerns the experimental situation here 0. We obtain very similar results 
for the currents within the BA. 

Let us compare these data with the experimental results. In Fig. 2 of Ref. |6[] one can 
identify the first resonant maximum very well before the domain formation sets in. It is 
located at U = 0.3 V, / = 0.64 mA. Dividing the voltage by the number of periods = 23 
we obtain eFd = 13 meV. Thus we have excellent agreement of both the position and the 
height of the peak. The low field resistance for the experimental I(V) curves is R ~ 570Q 
for the 35 well superlattice (Fig. 1 of Ref. JiT]) and R ~ 331Q for the 23 well superlattice 
(Fig. 2 of Ref. [0). Dividing by the number of wells we obtain d(eFd)/dI = 16Q and 
d{eFd)/dI = 14f2, respectively, which are in good agreement with our calculated value 
d(eFd)/dI = 12Q. The domain branches found experimentally exhibit maximum currents 
in the range between 0.35 mA and 0.55 mA which are lower than Im as common for field 
profiles where the high field domain located at the cathode. 

In Fig. we have also plotted some data points taken from Fig. 1 of Ref. |T7| (we divide 
the experimental voltage by the number of periods to obtain the effective field eFd). We can 
see that the onset of the second peak occurs for larger effective fields experimentally than 
in our calculation. This might be related to the fact, that a part of the voltage may drop 
outside the superlattice, so that the electric field in the sample is smaller than the estimation 
U/Nd used. Good agreement is found if we assume that a voltage drop U c = 0.33 V occurs 
in the contact. 



In Fig. 3 of Ref. Jl7j the experimental current is given in terms of e = 1.24 meV —U/N. 



The authors do not state that they found a maximum at U/N = 1.24 meV but just take 
this value from the measured intersubband spacing. Taking the voltage drop in the contact 
into account we have to use the value e = 124 meV—U c /N — eFd for our theoretical data. 
Then we obtain good agreement between the experimental data and the calculated curve as 
shown in Fig. [|. 

Analogously to Ref. 110] we calculated the current-voltage characteristic with domain 
formation using Eq. f| for the calculation of the current from one well to the next and the 
discretized Poisson equation. Here we assumed fast relaxation within the wells, so that the 
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upper subband is not occupied. Like in the experiment we consider a 23 well superlattice. 
For the electron density in well % we apply the boundary conditions n = = 0.95Nd/A. 
The resulting current- voltage characteristic for sweep up is shown in Fig. |5| which is in good 
agreement with the experimental findings. The high field domain is located at the anode in 
our simulation which is typical for the boundary conditions used and for sufficiently large 
dopings as can easily be understood within the general theory given in Ref. [[19] . 



IV. CONCLUSION 



We have developed a microscopic model to calculate the current in coupled quantum 
wells. The self-energies, which give the renormalized energies and widths of the spectral 
functions, are calculated directly from the nominal parameters describing the superlattice, 
without invoking additional adjustable parameters. The results for the current are in good 
quantitative agreement with the data from Ref. if we assume the barrier width to be 
1.5 nm less than the nominal value. We find good agreement for the low field behaviour 
including the position and height of the first maximum. Regarding the second peak the 
shape of the resonance is in good agreement if we assume, that 0.33 V of the total voltage 
drop outside the superlattice. Furthermore we have shown that the actual shape of the 
spectral functions due to the scattering seems to be responsible for a significant shift in the 
maximum of intersubband absorption. 
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FIG. 1. Square of the screened matrix element for impurity scattering. The full line gives 

2 2 

J2h Vh C (p) which determines the scattering in the BA. The dashed line gives Vq c (p) . The 
dotted line gives the same expression within the restriction to screening by electrons from the same 
well. The inset depicts the diagrams contained in the SSA. 
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FIG. 2. Spectral functions A a (Ek, E s ) for the first subband. The full and dotted line give the 
result within the SSA for E^ = and = 0.02 eV, respectively. The dashed-dotted and dashed 
line give the result within the BA for E^ = and E^ = 0.02 eV, respectively. The Energy E s is 
measured with respect to the bottom of the first subband. 
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FIG. 3. Current voltage characteristic (full line) for the parameter used in Ref. |17J except for 
b = 10.5 nm. The data points (crosses) are taken from Fig. 1 of [17]. Note that two different scales 
at the current axis are used here for the same curve. 
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FIG. 4. Current voltage characteristic near the second resonance. The crosses mark experi- 



mental data taken from Ref. [17]. 
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FIG. 5. Current voltage characteristic with domain formation for 23 wells for sweep-up of the 



voltage. The field distribution for U 
located at well number 0. 



1 V is shown in the inset where the injecting contact is 
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TABLES 



E a 


= 41.4 meV 


E b 


= 159.3 meV 


T a 


= -0.00539 meV (-0.0152 meV) 


R\ a 


= 1.61 • 10- 4 d (4.36 • 10~ 4 d) 


L P h 


= 0.38 meV 



TABLE I. Calculated parameters for the GaAs-Alo.3Gao.7As superlattices with well width 8 
nm and barrier width 12 nm. In brackets we have given the respective values for b = 10.5 nm. 
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